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Abstract. We address the open problem of existence of singularities for the complex Ginzburg- 
Landau equation. Using a combination of rigourous results and numerical computations, we describe 
a countable family of self-similar singularities. Our analysis includes the super-critical non-linear 
Schrodinger equation as a special case, and most of the described singularities are new even in that 
situation. We also consider the problem of stability of these singularities. 

1. Introduction. In this paper we study singular solutions to the initial value 
problem for the complex Ginzburg-Landau equation (CGL) 

du 

i— + (1 -ie)Au+ (1 +iS)\u\ 2a u = f, in IR'' x (0, T) , (1.1) 

u(x,0) = u (x) . (1.2) 

where u — u(x,t) is a complex- valued function denned in R d x (0,T) and satisfying 
suitable decay conditions as \x\ — > oo, the parameters e, S, a are non-negative real 
numbers, uq is a given initial condition, which is assumed to be smooth, with suitable 
decay as \x\ — > oo, and / = f{x,t) is a given function, also assumed to be smooth, 
with suitable decay as \x\ — > oo. 

We are mostly interested in the case e > 0, 5 > 0, and 2/d < a < 2/d + 1/2. 
However, some of our results are new even for the case e > 0, S > 0, a > 2/d, i.e., they 
also include the super-critical non-linear Schrodinger equation. Under the assumption 
e > 0, < a < 2/d + 1/2 it is possible to prove tha t fo r each smooth uq and /, with 
an appropriate decay at infinity, the problem ([Ti])-([l^) has a suitable weak solution, 
see [§] and |1^| . Such a solution is regular away from a closed set ScR^x (0, T) with 
V d ~ 2 / a (S) — 0, where V a denotes the parabolic a -dimensional Hausdorff measure, 
see [|l7|. It has been an open problem whether the singular set S can be non-empty. 

We present a very strong evidence, which is based on a combination of rigorous 
analysis and numerical computations, that singularities may indeed exist. We shall 
see that there appears to exist a countable set of different types of singular solu- 
tions. Among these solutions we identify (numerically) those which are stable under 
a suitable notion of stability defined in Section ^. It turns out that while most of the 
solutions are unstable, in certain cases there may exist more than one type of stable 
singularities. 

The case of the non- linear Schrodinger equation (NLS) (i.e., e = S = 0) and / = 
has been studied by numerous authors, see for example [|l(| || or the monograph |l^| . 
In particular, Zacharov (|p8|) conjectured the existence of self-similar singularities of 
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the form 

u(x, t) = (2/s(T - t)y^ +i ^ Q ((2k(T - i)r 1/2 M) , (1,3) 
where Q(£) is a complex valued function defined on (0, oo), with asymptotic behavior 

While no rigorous proof of this conjecture seems to be available, there is an over- 
whelming evidence based on numerical and formal analytical calculations supporting 
the existence of such singularities (see, e.g., jf| |?], Q). We also refer the reader to || 
for a rigorous result supporting the conjecture. 

In this paper we will argue that these singularities persist also for a certain range 
of e > and 5 > 0. In fact, we shall find many new self-similar singularities even for 
the case e = S = (and / = 0) . Using the self-similar singular solutions of the form 



( |l.3| ) one can easily construct singular solutions of fll.l] ) and ( [L2[ ) with compactly 
supported, smooth Uq and compactly supported, smooth /. 

From the form ([T^) of self-similar solutions one obtains the following boundary 
value problem for the function Q: 

(1 - ie)(Q" + ^—^Q') + i^Q' + i-Q - luQ + (1 + iS)\Q\ 2<7 Q = , (1.4) 

Q'(0) = 0, and (1.5) 
0(0 ~ C 1/a ~ luj/K as £ ^ oo . (1.6) 

In general, this problem does not have a solution for arbitrary values of parameters 



k, oj and hence the unknowns in (|1,4|)-(|1.6|) are Q, k and ui. 



When e = 0, see, for example, |Tj] for a discussion of fll.6|). For e > the condition 
Q(0 ~ ^-l/cr-iu)/K m £ n ^y j s dictated by the partial regularity result mentioned 



earlier, since u defined by ( |l,3| ) must be regular at almost all points of the form (x, T), 
x ^ 0. The presented results are based on a detailed analysis of the boundary value 
problem ( |l.4[) -(|L6|). The analysis employs both analytical and numerical techniques 
and is naturally divided into the following two steps: 

(i) First, we rigorously prove that for £i > 1 and for each set of parameter values e > 0, 

5 > and 2/d < a < 2/d+ 1/2 there exists a two-dimensional manifold of 
solutions to (|l.4|) on (£i, oo) with the correct asymptotic behavior at infinity. 

(ii) Second, we solve numerically the boundary value problem on (0,£i) with an 

appropriate (approximate) boundary condition for Q at the boundary point 
£ = £i. In this step the choice of £i is based on numerical evidence of 
convergence, and not on the rigorous estimates obtained in (i), which contain 
constants we did not try to evaluate exactly. 

We briefly summarize main results of our calculations. We note a (scaling) sym- 
metry in the problem Oj-(|uj): If A > 0, 8 e E and (Q(£),/t,w) is a solution of 
(L4)-(L : 6) then (A 1 / cr+l6 '(5(A£), A 2 k, \ 2 u) is also a solution. We chose a representative 



in each of these families of solutions by imposing suitable normalization conditions. 
We mostly work with the normalization u — 1 and ImQ(0) = 0, Q(0) > 0. Such 
solutions will be called normalized solutions and they are uniquely determined by two 
parameters (k, /z), where fi = Q(0). Some numerical results are better presented in a 
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different parameterization in which Q(0) = 1 is fixed, (k,lu) are used as parameters. 
We shall state the use of the latter normalization explicitly whenever it is used. 

We observe that for fixed values of d > 1, a > 2/d and e = 8 = there exists a 
coun table family of normalized solutions (Qj((,), Kj) — (Qj(Oi K j> J = 1 ; 2 , . . . of 
(|l.4[)-(p^). The j— th profile \Qj\, when extended to (—00,00) as an even function, 
has exactly j local maxima, and somewhat resembles the profile of the j— th state of an 
elementary quantum mechanical oscillator. A precise quantum mechanical interpre- 
tation of the solutions is more complicated, and is related to the so-called resonances 
or quasi- stationary states. The first solution of the family has been known, see or 
[Q, for example. We are not aware of any mentioning of the other solutions in the 
literature. All the solutions persist if e and 8 are perturbed to (small) strictly positive 
values. We now describe the behavior of these perturbed solutions for 5 — 0. We 
let Hj = Qj(0). As mentioned above, the solution (Qj,Kj) is determined by (Kj,fij), 
We observe that a branch of solutions parameterized by («j(e), Mi( e )) emanates from 



each point (Kj, jJhj). W e plo t the curves (e, Kj(e)), j = 1, . . . , 5, in Figure 3.1 for d = 1 



a = 2.3 and in Figure ^6| for d = 3, a = 1. In the other figures we plot the profiles 
|Qj(£)l of the corresponding solutions QJ(£) at certain points along each branch. 

An interesting feature observed in the behavior of the branches, is the existence of 
a turning point on each branch at e = e* . We see from the graphs that as e returns to 
zero along the branch, Kj(e) tends to zero, suggesting that the solution converges 
to a (radial) solution of the equation 

AQ-Q+\Q\ 2t7 Q = 0, mR d , (1.7) 

satisfying Q(x) — + as \x\ — > 00. 

Our computations presented in Section || clarify the structure of the diagram 
which turns out to be slightly more complicated than the picture suggested above. 
We conjecture the following: If d > 1, the solutions Qj corresponding to branch 
(ftj(e), Mj( e )) with an odd index j = 2k — 1 converge (as e and Kj(e) approach zero) 



to the k-th (normalized) radial solution of (1.7). We recall that the first of these 
solutions is usually called the ground state, and that for d — 1 there are no other 
solutions of dl.7| ) satisfying the appropriate boundary conditions. (See, for example, 
|Q for more details.) If d — 1 and also for j — 2k in the case d > 1, as e and Kj(e) 
approach zero, the profiles Qj separate into j approximate ground states which move 
away from each other. In particular, for j = 2k the profiles converge locally uniformly 
to zero. When d = 1 and j = 2k — 1, the profiles converge locally uniformly to the 
ground state. 



We tabulate results of our numerical calculations in Table |3,1| and Table |3.2 , 
We also looked at branches of solutions when 8 is related to e by 8 — re, with r > 
of order 10 _1 and 10°. The behavior was similar to the case 8 = 0, with the turning 
point e* getting closer to zero as r increased, as one might heuristically expect. 

Questions related to stability of the singularities are addressed in Section [|. Our 
calculations indicate that for the non-linear Schrodinger equation all the new singu- 
larities we found are unstable, and the singularity corresponding to («;i,/xi) is stable. 
The situation is more complicated for e > 0, see Section |] for details. 



Our interest in singular solutions to CGL stems from analogies between (LI) in 
the case d = 3, a = 1, e, 5 > and the three-dimensional Navier-Stokes equation 
(NSE). The two equations have the same scaling properties and the same energy 
identity. Moreover, the existence and partial regularity theory of weak solutions for 
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NSE and CGL are similar (with CGL being technically easier), see |17j]. The analogy 
between NSE and CGL may be rather superficial and may break down at any deeper 
level. However, at the same time there are no known properties of solutions to the 
Navier-Stokes equation which would prevent the same scenario as presented here for 
CGL. 



The formula for the Navier-Stokes equation corresponding to (1.3) would be 



u(x, t) ~ (2k(T - t))- 1 ' 2 U ((2k(T - i))- 1/2 x, t) , (1.8) 

where r = j- In , and U is a suitable divergence-free vector field periodic in r with 
suitable decay in the self-similar variable y = (2n(T — £)) -1 / 2 x. 

The case, dU/dr = 0, was already considered by Leray It was proved in |llj] 
and i n gre ater generality in [[jj) that NSE does not admit non-trivial solutions of the 



form fll.q ) with U independent of t. The problem is open for U periodic in r. 

We finish the introduction with the following speculation. Most of the singulari- 
ties we have found are unstable, hence it is unlikely they would be observed in direct 



numerical simulations of the initial value problem (|l.l[)-(1.2) or in physical experi 



ments that are modeled by CGL. Could it perhaps be the case that NSE does admit 
singular solutions (say of the form ( |l,§l )), but all of them are unstable and therefore 
more or less impossible to be detected in direct numerical simulations or physical 
experiments ? 

This intriguing scenario was once suggested to one of the authors by Sergiu 
Klainerman during a lunch-break conversation at a conference in Southern California. 

2. Analysis of the profile equation at infinity. In this section we study 
solutions of (1.4) in the interval (£i,oo) satisfying the condition Q(£) ~ £- 1 / a -' lu, / K 



as £ — * oo. Heuristically one expe cts that the behavior of such solutions is mainly 
governed by the linear part of (1.4): 



(1 - ie)u" + (1 - ie)—- —u' + in^u' + —u - uju = (2.1) 



Equation (|2.l|) is equivalent to Rummer's equation, also known as the confluent hy- 
pergeometric equation. The solutions of this equation are well-understood, see for 
example Jl3), and one can hence get a more or less complete picture of the behav- 
ior of solutions to ([Tl]) . We shall use analytical tools from the theory of confluent 



hyper-geometric equations to describe solutions of the full equation ( |l.4|) . 
A canonical form of Rummer's equation is 

d 2 w x dw 

z— r + ft-z- aw = 0, (2.2 

dz z dz 

and the equation ( |2.l[ ) is transformed into this form by letting 

—in £ 2 1/1 iu\ d , 

There is voluminous literature on this equation and properties of special functions 
(confluent hyper-geometric functions) which appear as its solutions. We recall some 
of the properties and, for the convenience of the reader, we also sketch how to derive 
them. For more details about confluent hyper-geometric functions we refer the reader 
to @, §, @. 
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A classical formula for a solution of (2.2) is given by 



1 f°° 

U(a,b,z) = —- e- u ^- 1 (l + t) h - a - 1 dt, (2.4) 
1 \ a ) Jo 

where T is Euler's gamma function. The integral is clearly well-defined for Re a > 
and Re 2 > 0. The factor 1/T(a) is not essential for our analysis, nevertheless, we 
include it to keep our notation in agreement with the standard one. The role of this 
factor is to normalize the leading term in the asymptotic U(a, b, z) = z~ a (l + 0(z -1 )) 



as z — > oo. It is easy to check by direct calculation that the function U given by (2.4) 



solves the equation (2.2). We have 



£ j -U(a.,b : z) = f ^- / (-/>''< '■■/"■ '(1 +/)''-"-'<//. 



and after substituting to (2.2) we see from a simple integration by parts that the 
equation is satisfied. 



By letting zt = s in (2.4) we obtain (for Re z > 0, Re a > 0) 



If 00 / s \6-a-l 

U(a,b,z) = ——z~ a e~ s s a - x (1 + -J ds . (2.5) 
r(a) J \ zJ 

The above expression is used to extend the definition of U for Re a > and z € C with 
— 7r < argz < 7r, since the integral is convergent and an analytic function of z under 
these assumptions. By formally expanding the term (l+s/z) f, ~ a_1 and integrating the 
resulting (formal) series term-by-term, we obtain the following asymptotic expansion 

U(a, b, z) = (j2 (a)fc(1 ± a ~ b)k (->)-* + O(N-)) , (2.6) 

where (a)o = 1, (a)fe = a(a + 1) . . . (a + fc — 1). The formal calculations of the 
asymptotic expansion ( |2.^) can be easily justified rigorously by splitting the integral 
in ( |2.5| ) as J °° = j^ 1 + j g with Si = \z\ 1 ^ 2 , for example. 

One can check easily by a direct calculation that the function 

V(a,b,z) = e z U(b-a,b, -z) (2.7) 

is another solution of Rummer's equation and that the functions U and V are linearly 
independent. The Wronskian of U and V is given by 

uW_ v w =e±iv{b _ a)z _ be ^ 

dz dz 



is 



where the sign + is for Imz > and — in the opposite case. (Formula (^_ 
easily derived from the fact that the Wronskian satisfies the differential equation 
y' + (b/z — l)y — and from the asymptotic expansion (0).) 

We need to know the behavior of U and V in the region — n/2 < argz < 0. 
Taking into account the definition of V, we see that it is sufficient to control U in 



the region — n/2 < argz < n. Formula (2J5) is suitable for analysis in this region if 



Re (b — a) > since the integral in ( |2.5j ) is the uniformly absolutely convergent when- 
ever z approaches a point in (— oo, 0) from the upper half-pl ane. In our applications 
the condition Re (b — a) > is always satisfied and therefore ( |2.5j ) is sufficient for our 
analysis. 
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We now have sufficient information about the solutions of (2.2), and hence also 



( p.lD , to be able to proceed with the analysis of the inhomogeneous equation 

. s d- 1 , 



{1 - ie)u" + (1 - ie)- 



- a + in^u H u — oju = f(£) 

a 



(2.9) 



-1 / a—iu) I k 



as 



for £ G (£i, oc). We assume that the function / is decaying sufficiently fast as £ — > oo. 
We are interested in solutions of (2J3) which have the asymptotics u ~ £ 
£ — > oo. We denote P, £7 two linearly independent solutions of (1.4) 

P(e) = P(K,u>,e;t) = U(a,b,z), 



where z = 



-in e 

1 ^ie~2 



2 \cr 



b = 



The parameters d and tr are held fixed in the perturbation analysis described below, 
therefore we do not indicate the depende nce of P and E on them. The Wronskian 
W = PE' - P'E is easily computed from 



W(e)=W(o;,«,e;0 



= ±7ri(b— a' 



1 -ie 



(2.10) 



The next step is to use the standard variation of constant to obtain solutions to 
(2.9) satisfying u(£) ~ ( e- 1 / cr -W K as £ _> oo. We write the solution in the form 



u(0=c 1 (£)P(0+c 2 (0E(0, with 
« / (0=d(0i >/ (0+ca(0^ , (0. 

and C2(£) ^ "sufficiently fast" as £ — > oo. We obtain 



(2.11) 
(2.12) 



(2-13) 



which together with the condition C2(£) — > at the infinity gives a formal expression 
for the solution 



«(0 = 7^(0 
where 7 6 C is a constant and 



(2.14) 



-^^(O^tl)^"^) for £1 < 77 < £ 
-{irsj^O^WW-^T?) for£<77 



The strategy for finding solutions of the non- linear problem (1.4) in (£1, 00) with the 
required decay w(£) ~ £-i/<t-iu/k m £ n ^y should now be clear. Heuristically we ex- 
pect that the manifold of such solutions will be a deformation of the one-dimensional 
complex subspace { 7 P | 7 £ C}, at least in a neighborhood of the origin. The defor- 
mation is from 7P to the fixed point of the operator T 



Tu(0 = 7 P(0 



(l + i6)K(£,r 1 )\u(ri)\ 2 ' T u{T 1 )dr). 
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The outlined strategy can be successfully carried out by using the properties of 
Kummer's functions recalled above. The fixed point theorem can be applied in the 
Banach space 

X* = {u G C([£i,oo)) | sup |e| 1/ff -*|u(OI < 00} 



equipped with the norm 



H, = sup l^l^KOI. 



110 

€>6 



This approach is obviously standard. However, there are some subtle points in the 
situation studied here due to oscillatory behavior of the function E, see the Appendix. 

The main result of this section, which can be derived in a fully rigorous way from 
the above analysis is the following theorem. (A complete proof of the theorem is 
presented in the Appendix.) 

Theorem 2.1. Assume 1 < d < 3, 2/d < a < 2/d + 1/2, < m < k 2 , w\ < lu 2 , 
< ei, Si < S 2 . 

There exists pi > such that for each £i > 1 and each 

((3,K,w,e,6) G A P1 = {/?GC| \/3\ < pi} x [ki,k 2 ] x [w 1 ,w 2 ] x [0,ei] x [6 1 ,5 2 ] 
the boundary value problem 

(1 - ie){Q" + ^—^Q') + inZQ' + i-Q - wQ + (1 + iS)\Q\ 2a Q = , (2.15) 

Q(£i) = /3, «nd (2.16) 
Q(0 ~ C 1/a ~ iu/K as £ -» oo . (2.17) 

considered in (£i,oo) /ias a solution 

F(0=F(/3,K,u,e,6,ti;t). 

Moreover, F can be constructed in such a way that the following conditions are satisfied 

(i) The mapping from A Pl to X$ defined by 

(/?, k, uj, e, 5) i > k, w, e, 5, £ i; •) 

is C up to the boundary for each d > 0. 

(ii) The complex-valued function defined by 

dF 

(/3, k, u, e, S) i-> ^ (/?, «, w, e, (5, £i, £i) 

is C 1 (Mp to the boundary) in A P1 . 
(hi) F and its derivatives F^> have the following asymptotic expansions 

F(£) = (jlaiC 21 + 0(r 2(rl+1) )^) , (2.18) 

F w#) = &L U- 1/a ~ iw/K it ai Z) + o(r 1/CT - 2( " +1) - fe ) • (2.19) 
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(iv) We have F(0, k, uj, e, S, £1, £) = and 



— (0,/t,w,e,5,£i,£) - 



9/5-' — — - P(«,w,e,Ci) 
Proof: See the Appendix 



Remark 2.1. For e > i/ie function u(x, t) given by has to be regular at all 
points (x, t) wit h x ^ by the partial regularity theorem proved in fiT\] . Therefore any 



solution of nLl )-(L6) must admit an asymptotic expansion of the form stated in (Hi) 
of Theorem We note that the convergence of the series a i£,~ 21 * n the asymptotic 
expansion of Q is equivalent to the analycity of u (in t) at the points (x,T), x ^ 0. 
The asymptotic expansion of the profiles Q does not converge and therefore u is not 
analytic in t at any point (x,T). 



Theorem 2.1 together with elementary perturbation arguments can be used to 



show that if some non-degeneracy conditions are satisfied, then every solution of 



(1.4)- (1.6) for e = 5 = will persist (with a slight deformation) for small e > and 
S > 0. We will briefly describe this standard procedure for the convenience of the 
reader. 

First we consider solutions on the finite interval (0, £i] to the initial value problem 

(1 - ie)(Q" + ^—^Q') + w£Q' + i-Q - loQ + (1 + iS)\Q\ 2 °Q = , (2.20) 
Q(0) = m, Q'(0) = 0. (2.21) 

We denote the solution (if it exists) by G(£) = G(fi,K,uj,e,5;£). Clearly the set 
of parameters (/U, K, uj, e, 6) for which G is well-defined is open. We define 

k, uj, e, S) = G(ji, K, uj, e, 5; £i) . (2.22) 



Assume that the boundary- value problem (|l.4|)-(|l.6D has a solution which sat- 
isfies Q(£i)= [i. With a slight abuse of notation, let us denote such a solution by 
Q(p, k, uj, e, 5; £). Clearly Q(/i, k, uj, e, 5; £) is defined only on a submanifold of the 
parameter space, but this will not be important in what follows. As we have seen in 
the introduction, we have 

A^+^QOx, K, uj, e, 6, A£) - Q(A 1 / ff+l V, A 2 k, X 2 lu, e, 6; £) , 

for all A > and 8 £ [0, 2ir). Therefore we can work with normalized solutions, i. e. 
we assume that uj — 1 and that Q(0) is real and non-negative. Assume > 0, n > 
and suppose Q(fJ-o, K o, 1, 0, 0; £) exists. 
We set 

OG dF 
g(fi, k, e, S) = «, 1, e, 6; £i) - ^ (Z 3 ^ K > 1 : e > S )> K > 1 : e > S > £i> £i) ■ ( 2 - 23 ) 



where /3(^, n, uj, e, S) is defined by fl2.22|) . By our assumptions and by Theorem the 
mapping g is well defined and continuously differentiable in a set of the form [pi,^) x 
(uji,ui2) x [0, t\) x (8\,52) containing the point (/io, Kq, 0, 0). Since g(po, no, 0, 0) = 0, 
we see that the equation 



g(fx,K,e,6) = 
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has solutions for small e > 0, 5 > if the following non-degeneracy condition is 
satisfied: 



det 



4 H j =M) at. (//„. /.-„.(). (I) 



(2.24) 



where g 1 — Reg, g 2 — Img and subscripts denote partial derivatives with respect to 
the corresponding variables. 

Thus a non-trivial solution of ( |l,4| )-( fr6| ) for e = 5 = also gives a solution of 
dl.4|)-(p^) for e, 5 > if (2.24) is satisfied. Based on our numerical calculations 



described in the next section, we conjecture that (2.24) is satisfied for a > 2/d. 

3. Numerical results. Theorem 2.1 allows us to rewrite the boundary value 



problem (1,4)-(1.£) as a boundary value problem on a finite interval (0, £i) in the 
following way. 



(1 - ie)(Q" + ^—^Q') + ik£Q' + i-Q-Q + (l + iS)\Q\ 2a Q = . 

Q'(0) = 0, 

<?(&) = 0, 
dF 

Q'(a) = 7^03,K,W,6;£i), 



(3-1) 

(3.2) 
(3.3) 

(3.4) 



where the unknown quantities are Q, (3 and k. Of course, the problem (|3.l|)-(|3.4|) 
is equivalent to the equation g = in the previous section. If w e a p prox imate F 
by the first few terms of its asymptotic expansion, the problem (3.1)-(3.4) can be 
solved numerically. In our numerical computations we investigated the dependence 
on £i and on the number of terms of the asymptotic expansion of F. It turned out 
that £i ~ 30 and the first term of the asymptotic expansion already worked very 
well. However, many of our computations were done with the first two terms of the 
asymptotic expansion (2.18). In the case e = 8 = Q, i.e. NLS, the value of (fit, k±) was 
computed in j?J using completely different approach. The values presented in that 
paper are in an excellent agreement with our computations, see below. 

The first term of the asymptotic expansion of F is F ~ /?(£/<!;i)~ 1 / ' -,w / K . Using 
this approximation, we obtain from (EO)-(3.4) 



(1 - ie)(Q" + ^—^Q') + iKZQ' + i-Q - 
Q'(0) = 0, 
6Q'(6) 



= 0, 



a k 



(3.5) 
(3.6) 
(3.7) 



Higher order approximations can be derived in a similar way. Note that in the formu- 
lation of (3.5) we already fixed the normalization oj = 1, so that the unknowns are Q 
and k. The boundary condition (3.7) is also closely related to the boundary condition 
used in Hi, || for simulations based on solving time dependent problem in the PDE 



(yjH 



There are essentially two approaches to the numerical solution of the boundary- 
value problem (3.5)-(3.7). One can use collocation methods to approximate the 
boundary-value problem and then to apply Newton's method to the discretization 
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of the non- linear operator which defines the equation (3J5). Implementation of this 
strategy requires further changes in the formulation since the linearization of the non- 
linear operator always has zero in its spectrum due to the S' 1 -equivariance of the 
equation. Therefore a shooting method was easier to implement and it also proved 
to be sufficiently accurate. Because of well-known sensitivity of shooting methods 
to problem parameters we performed computations in different normalizations: with 
Q(0) = 1 fixed and parameters (k,uj) as the unknowns as well as with w = 1 fixed 
and parameters (k, /i) as the unknowns. Moreover we compared both backward and 
forward shooting methods on the interval (0,£i). All computed solutions turned out 
to be in a very good agreement. 

Here we describe only the shooting method using the normalization w = 1, in 
which we calculate + (l/a + i/n)Q(£i) as a function of the parameters fx = 

Q(0), K, e, and 8, i.e., we solve the equation 

f{H, k, e, 5) = 6Q'(6) + (1/* + %/k) Q(6) = . (3.8) 

The integration of the underlying ODE must be done with sufficient accuracy. We 
compared various ODE solvers. A variable-order, variable-step Adams method as 
implemented, for example, in the NAG library proved to be sufficiently accurate in 
most of the calculations. To locate the initial values for each branch we inspected 
the two dimensional subspace of the parameter space given by (/j,, K, 0, 0) for d = 1 
and cr = 2.3 and computed the degree of the function / restricted to that subspace 
along various curves. The other solutions were then calculated by continuation. This 
was done with the help of bifurcation analysis package developed as a part of [ |l2| . 
The implementation of the path-following procedure with a Newton corrector step 
can be done efficiently as the linearization along a solution is evaluated at the same 
integration step as the solution. 

As we described in the introduction, it appears that the equation f(fx, K, 0, 0) has 
countably many solutions (Kj,fij) in the region k > 0, fi > 0. The corresponding pro- 
files \Qj\, when extended to (— oo, oo) as even functions, have exactly j local maxima. 

Remark 3.1. As mentioned in the beginning of this section the solutions (Kj,fij) 
may slightly depend on the value £i and the approximation of F. For solutions de- 
scribed here we tested the dependence of the results on £i for £i € [20, 100], and also 
on the approximation of F by taking either one or two terms in the asymptotics ex- 
pansion of F. We also directly compared our numerical solutions in intervals of the 
form (t)i,£i) with the explicit formulae given by one or two terms of the asymptotic 
expansion of the solution in (?7i,£i) for various values 0/771. All these tests indicated 
a good convergence of our approximations. Based on these tests, we estimate that the 
error in the values of the "roots" (kj, fXj) is of the order 10 -3 or better. 
Case I (d = 1, a = 2.3, 5 = 0): Results for the one-dimensional case are tabulated 



in Table 3.1 



The continuation of solutions parameterized by e is depicted in Figure gjj. To 
give the reader a good idea about the form of the solutions we plot profiles \Q(£,)\ a t 
a few points on each branch (at the point e = (solution to NLS), at a point on the 
upper part of the branch and at another point on the lower part of the bran ch) . We 
used both normalizations, Q(0) — 1 in Figure 3^-3^, and u) = 1 in Figure 3.4 - 3.1! 



Normalization Q(0) = 1 is convenient for computing solutions along the odd branches 
as these solutions exhibit a maximum at the origin which is also present in the solution 
for £^0on the lower part of the branch. 
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12388 





20263 
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26678 
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47127 
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76651 


5 





03455 





21643 


1 


07947 





15225 





70345 


6 





03099 





18185 





92714 





25750 


1 


41624 
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02803 





15667 


1 


05430 





12284 





78409 
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02559 





13756 





95061 





17365 


1 


26236 



Table 3.1 

Case I (d = 1, a = 2.3, <5 = 0): e* defines the position of the turning point; coordinates (k, p.), 
where fi = Q(0) define the initial point of each branch, when e = 0, in the normalization (a = 1. 
The coordinates refer to the same solutions in the normalization Q(0) = 1. 



"ir3 




Fig. 3.1. Case I (d = 1, a = 2.3, 8 = 0): solution branches (e, Kj(e)); normalization Q(0) = 1 
(left) and u> = 1 (right). Turning points are denoted by cu The solid line indicates stable solutions, 
the dashed line indicates unstable solutions. See Section y| for details concerning stability issues. 



We now describe the behavior of solutions along the branches. We observe that 
as e returns to zero after passing through the turning point, we have — » 0. 

For j = l the solution Qf(£) approaches a specific solution, usually called the 
ground state, of the equation 

-u" + u=\u\ 2a u. (3.9) 

For j > 1 the profile |Qj(£)| seems to separate into j copies of the ground-state 
solution which are moving to infinity. For j even all of them "escape" to infinity, while 
for j odd one will stay at the origin and the rest will move to the infinity as e — > 
and Kj(e) — > 0. 

Case II (d = 3, a = 1, S = 0): Results for the three-dimensional case are tabulated 
in Table |J. 

Some data for NLS (e = 6 = 0) and the "basic solution" (corresponding to the 
beginning of our first branch) are available in the literature and can be used to estimate 
accuracy of our calculations. One can see a very good agreement of our solution for 
j = 1 with values n — 0.917 and = Q{0) = 1.885 obtained in from simulations 
that used the dynamical rescaling method applied to the initial value problem (1.1)- 
& 
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Fig. 3.3. Case I (d = 1, a = 2.3, <5 = 0): profiles |Q(£)| at points along the j-th branch; j = 5 
fl e ft)> j = 7 (right); normalization Q(0) = 1. 



The behavior of solutions when e returns to zero is qualitatively similar to the 
one-dimensional case, but there are also some new interesting features. 

First, the maximal possible e* is not attained on the first branch (j — 1). Second, 
the behavior as e — > and Kj(e) — * is different for even and odd values j. For 
j = 2k — 1 (odd) the solutions QJ(£) converge (along the jth-branch) to the fc-th 
radial solution of the problem 

-Au + u = \u\ 2a u in R 3 , (3.10) 
w^O as \x\ — > oo . (3-11) 

For j = 2k (even), one observes a similar behavior as in the one-dimensional case. 

Remark 3.2. From the regularity theory for CGL (see fify) one expects that the 
value of e cannot cross a certain finite value e, as we follow solutions along any branch. 



What we see in the calculations is a turning point on each branch. In Table 3.1 and 



Table 3.i we list for each of the calculated branches the parameter e* which is the 



maximal value of the perturbation parameter e reached along the branch. 
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Fig. 3.5. Case I (d = 1, <r = 2.3, 5 = 0): profiles |Q(£)| a * points along the j-th branch; j = 3 
(left), j ' = 4 (right); normalization u) = 1. 



4. Stability of singular solutions. A natural question regarding the stability 
of the self-similar singularities constructed in the previous section is for example the 



following: What will be the behavior of solutions to (|1.1|) and (1.2) if uq is a slight 
perturbation of a solution Q of the profile equation ( |l.4|) -(1.6) and / is small? (One 
can think for example of taking uq = ipQ + ut, where <p is a compactly supported 
smooth cut- off fu nct ion which is identically one in a large ball, Q is a non-trivial 
solution of (|T|)-(|l|), u\ is small and compactly supported, and at the same time 
taking / which is small and compactly supported in x.) 

Here we present an approach to this problem which uses the method of dynamical 
rescaling (see for example Jl(J), which seems to be natural in this context. We empha- 
size that we will not obtain fully rigorous analytical results which would completely 
answer the question raised above. Our goal is to present some preliminary calcula- 
tions which seem to be adequate to the issue of stability in the context of numerical 
simulations. 



We briefly recall the main idea of the method of dynamical rescaling. We consider 
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93792 


5 


0.18168 
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Table 3.2 



Case II (d = 3, a = 1, S = 0): e* defines the position of the turning point; coordinates (k,/i), 
where fi = Q(0) define the initial point of each branch, when e = 0, in the normalization u = 1. 
The coordinates (k,u>) refer to the same solutions in the normalization Q(0) = 1. In Km the solution 
corresponding to j = 1 and e = was calculated by the method of dynamical rescaling, which gave 
k = 0.917 and fi = 1.885. 




Fig. 3.6. Case II (d = 3, <r = 1, <5 = 0): solution branches (e, Kj(e)); normalization Q(0) = 1 
(left) and ui = 1 (right). Turning points are denoted denoted by o.The solid line indicates stable 
solutions, the dashed line indicates unstable solutions. See Section U for details regarding stability 
issues. 



a solution u : M. d x (0, T) — > C of the complex Ginzburg-Landau equation 

i— + (1 - ie)Au + (1 + iS)\u\ 2a u = . (4.1) 

Motivated by the scaling invariance of solutions 

u(x,t) -> X 1/a u{Xx,X 2 t) 
of d4.l|) , we write the solution u as 

u(x,t) = L- 1/a (t)v (L- 1 ^^) , 
where L(t) > is to be chosen and rfr = L~ 2 {t) dt. From ( fi~l| ) we obtain 

i^- + «(r) + -v) ie)Av + (1 + i5)\v\ 2 °v = , (4.2) 



dr \ 9^ a 

where k(t) = -L _1 (r) dL(r)/ dr = -L(t) dL(t) / dt. 

We now chose L(t) so that, roughly speaking, the typical length-scale over which 
v oscillates is 1. If u develops a singularity at t — T (and it is regular in (0,Ti) for 
each Ti < T), then L(t) > in (0, T) and L(t) -> as t T. 
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Fig. 3.8. Case II (d = 3, a = 1, 6 = 0): the profile |Q(£)| at points along the j-th branch, 
j = 1 (left), j = 2 (right); normalization oj = 1. 



One way to control oscillations of v is to impose a condition J{v{., r)) = 1, where 
J is a suitable functional controlling regularity of solutions to ( |4.2D . One is then led 
to the following system 

+ w(t) + -v) + (1 - ie)Av + (1 + i5)\v\ 2 °v = , (4.3) 



dr \ d£ a 

J(v(.,t)) = 1, (4.4) 

where the unknowns are v and k(t). Various choices of the functional J have been 
used for numerical calculations. For example, in [jio], || the functional 

Jo(v) = [ \Vv\ 2 dx 

was used. Although it does not directly f ollow from the known regularity theory that 



J7q controls the regularity of solutions to (4.2) in the case e = 5 = this choice turned 



out to work satisfactorily in the numerical computations. 

In the following analysis we assume that the functional J is invariant under the 
action of the symmetry group S 1 , i.e., J{e ie v) — J{v) for all 9 £ [0, 2ir). Self- 
similar singularities of the form ( |l.3| ) correspond to solutions of ( |4.3| ) of the form 
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v(£,t) = e tor Q(^), i.e., to S^-orbits of solutions to the problem (0)-0)- We will 
consider the linearized stability of these orbits, which seems to be the simplest natural 
notion of stability in our context. Suppose (i>o(£,t), kq(t)) is a solution of (4J3) such 
that wo(C> r ) — eV *° T Q(Q and kq(t) = kq = const. We consider a perturbed solution 
t), k(t)) in the form 

v&T)=e i "° T (Q(0+w(Z,T)) 
k(t) = k + x{t) , 

where w and >c are infinitesimally small. A simple calculation gives 

'''- =Cw + k{t) (zQ' + -q) , (4.5) 



dr ' ' V a 

{J'{Q),w) = Q, (4.6) 



where the operator C is defined as the linearization of the equation (4.2) along the 
solution (i>o,ko), i- e -' 

Cw = (i + e)Aw - k + -w) - iuj w + i(l + a)\Q\ 2a w + ia\Q\ 2,T ~ 2 Q 2 w 



and J' denotes the derivative of J . (Note that the operator C is not complex linear 
and therefore it is natural to carry out the analysis in the real representation.) 

There are two eigenfunctions of C that can be formally derived directly from the 
invariance of the profile equation ( |l.4| ) under the scaling symmetries 

(Q(0,K,w)-*(e"Q(0,K,w), (4.7) 
(0(0, «, to) -> (A^+^Q^), A 2 k, A 2 c) . (4.8) 



The invariance under the S^-symmetry (4.7) implies existence of an eigenfunction 
Y\ = iQ with the eigenvalue zero, while the symmetry under (4.8) leads (formally) to 
an eigenfunction Y2 = £Q' + (l/a + iujo/no)Q with the eigenvalue 2ko- We note that 
Y 2 {£,) = O^- 1 ^- 2 ) as £ -> 00, due to ( |2~18| ) and ( pj|) . We have 

(4.9) 
(4.10) 



CY X = , 
£F 2 = 2k F 2 
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We rewrite (4.5) as 



OT \ Kq 

(J'(Q),w} = 0, 



(4.11) 
(4.12) 



where the unknowns are w(£,t) and x(r). 

Proceeding further with our formal reasoning, we view ( 4.11 ) as a dynamical 
system with one linear constraint in a suitable linear space W of functions on M. d . In 
this formal analysis we will not try to specify W. For a rigorous analysis it would 
be natural to try to find a suitable Banach space containing all smooth compactly 
supported functions, together with the functions Y\ and Y2. Suppose that the space 
W can be decomposed into a direct sum W = MYi®WY2®Z where Z C Wis invariant 
under C. In the finite-dimensional situation, a sufficient condition for this would be 
that the eigenvalues and 2kq are simple. Assuming that such a decomposition exists, 
we write 



w(£,t) = wi(t)Y 1 (Z)+w 2 (t)Y 2 (0 + Z(£,t) 

where w\, w 2 are scala r func tions and Z (.,t) £ Z. 

After substituting (|4.13|) into ([4.1l|) and using (J'(Q),Yi) = 
the invariance of J) , we obtain 



dz 

dw 2 

W 2 



-x{t) — 
2n w 2 + x{t) 

(J'(Q),Y 2 ) 



(4.13) 

(a consequence of 

(4.14) 
(4.15) 
(4.16) 
(4.17) 



where we assume that (J'{ Q), Y%) ^ 0. T his is a completely natural assumption in 
the context of (4.2). From (4.14) - (4.17) one can see that in a finite-dimensional 



situation and under the assumptions stated above, the following conditions would be 
equivalent: 



(i) If (w,>c) is a solution of (4.5), then u>(£,t) approaches exponentially aYi(£) for 

some a £ K and >f(r) approaches exponentially zero. 

(ii) The spectrum of C\z belongs to the set {z £ C | Re z < 0}. 

In a finite-dimensional situation and under the assumption that the eigenvalues 
and 2kq are simple, condition (ii) would be equivalent to the condition that all the 
spectrum of C except and 2kq lies in {z G C | Re z < 0}. 

The condition (i) is exactly the linearized orbital stability of the orbit (e luJ ° T Q, Kq) 
for the system (|0|), (jDj). 

It seems to be a non-trivial problem to put the above formal analysis on a rig- 
orous basis in the infinite-dimensional setting. However, the formal analysis strongly 
suggests tha t in a finite dimensional situation which arises in numerical approxima- 
tions of J4.3| ) the stability of the solutions (e lu ° T Q(£,), kq) should be governed by the 
spectrum of an appropriate approximation of the operator C. 

We remark that our stability analysis is independent of J, except for the natural 
assumption (J'(Q),Y 2 ) ^ 0. 
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We numerically calculated the approximation of eigenvalues for a discrete ap- 
proximation of C in the space of radial fun ctions along the branches of solutions 
parameterized by e. The original system ( |4.3| ) was truncated to a finite interval (0, £1) 
by imposing the time-dependent boundary condition at £ = £1 

^Ki.^ + ^l^r) + ^v(£i,t) = 0, (4.18) 
which is a time-dependent equivalent of the bo unda ry condition derived from the 



asymptotic expansion in Section g. The condition fl4.1q ) was also used in the numerical 
simulations in |Q| and ||. The computations were carried out for £x = 30. 

The accuracy of our numerical approximation can be checked indirectly by com- 
paring the predicted eigenvalues and 2kq with the corresponding eigenvalues we 
obtained from our numerical calculations. We saw a very good agreement, in most 
cases the error was of the order 10 -4 . Both and 2ko appeared simple, except in 
some natural degenerate cases when other eigenvalues were crossing them as we moved 
along branches. 

Case I {d = 1, a = 2.3, 5 = 0): The calculations confirmed what one intuitively 
expects: The solutions on the upp er p art of the branch j = 1 are stable and all other 
solutions are unstable. In Figure |3.l| we used a solid line for the stable parts of the 
curves and a dashed line for the unstable parts. 

Case II (d = 3, a = 1, S = 0): The situation is similar with one notable exception. 
In our computations we detected stable solutions also on the lower part of the branch 
j = 2. Accuracy of our approximation did not allow us to decide whether all solutions 
on the lower part of the branch j — 2 for < e < e* are stable, since as e approached 
zero we observed some eigenvalues very close to the imaginary axis. In Figure 3.6 we 
used a solid line for the stable parts of the curves and a dashed line for the unstable 
parts. 

We emphasize again that the calculations we carried out only deal with stability 
in the space of radial functions. Based on computations in ||, it appears that the 
solutions on the upper part of the branch j = 1 are also stable with respect to 
perturbations that break the radial symmetry. 

It is clear that further work is required to fully clarify issues concerning the 
stability of solutions described in this paper. 



2.1 



5. Appendix. In this section we give a full proof of Theorem 
We fix r% > and $ > 0. Values of these parameters will be chosen later. We 

denote by X$ the Banach space of continuous functions u : oo) — » C for which the 

norm 

M* = sup \£,\ 1/a ^\u(0\ 

is finite. In this proof we use the following notation 

A = [«i,«2] x [wi,w 2 ] x [0,ei] x [Si,5 2 ] 
A = {7 e C| |7| < ri} x A. 

A point in A is denoted by A = (Ai, A 2 , A 3 , A 4 ) = (k, ui, e, S) and similarly A G A rep- 
resents (Aq, Ai, A2, A3, A4) = (7, k, ui, e, S). The functions P, E, W, and K introduced 
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in Section | are written as P = P(X, £), E = E(X, £), W = W(X, £, rj), K = K(X, £). 
With a slight abuse of notation we sometimes write also P = P(A, £) even if the 
function does not depend on Xq. 

Lemma 5.1. There exists C > such that for A e A, £, ?y > £o we have 



\K(X^,V)\ < 



C £-d+l/* v -l-l/<r+d f Qr Z <£< v . 



(5.1) 



Proof: This statement is an easy consequence of the definition of K and of the 
asymptotic expansion for U a s z — > oo in the sector — tt/2 < argz < it. An important 
point in connection with (2.6) is that the constant in the remainder 0(|z|~") can be 
taken same when the parameter a runs through a compact subset of {z | Rez > 0}. 
(This can be seen for example from (pH|).)D 

Lemma 5.2. Assume (2a + l)i9 < 2 + 2/cr - d. Then the formula 



T(X,u)(0=^P(X,0 



(l + iS)KCX,^v)Hv)\ 2a u(r,)df 1 



Si 



defines a continuous mapping T : A X X& — ► X& . 
Moreover, there exists C > such that 

|T(A,u)|,<C| 7 |£r*+CS 



2+2aS 



2a+l 



(5.2) 



(5.3) 



\T(X,u) T(\v)l < C^ 2+2 °»\\u - vl (\\u\\f + \\v\\l c 



(5.4) 



for all X 6 A and all u, v G X# . 

PROOF: The convergence of the integrals in (5.2) and ( |5.3[ ) follows easily from 
Lemma 5.1. To get the estimate (5.4) we use Lemma 5.1 together with the elementary 
inequality 

|N 2 ^i - M 2 ^ 2 | < M\zi - z 2 \ (|z x | 2ct + |z 2 | 2ff ) , 

which holds for all z%, z% € C if M > is properly chosen. 

The continuity of the mapping T can be proved by suitably splitting the integral 
over (£i,oo) into two integrals, one over and the other one over (£2,00). We 

can then estimate the integral over (£2, 00) by using Lemma 5.1. □ 



The estimates fl5.3| ) and (5.4) show that T(A, .) is a contraction of the ball B p 
{u 6 Xg I \u\$ < p} into itself if p is such that 



Cr^ + C^ 2+2 ^p 2a+1 <p, 
2Cp 2 °^ 2+2 ^ < 1 . 



(5.5) 
(5.6) 



Proposition 5.3. Assume (2er+l)tf < 2+2/ a-d. If ([5.3j) and t \5.(\ ) are satisfied, 
then for each A £ A the mapping T(X,.) has a unique fixed point u(X) — u(A,£) in 

Bp = {u e X# I \u\# < p}. 

Moreover, the mapping X 1— > u(X) is a continuous map ping from A to X$. 
PROOF: The proposition follows directly from Lemma I5T2] and the Banach Fixed 
Point Theorem. □ 

Next we need to establish that the solutions u(X, £) have the required regularity 
properties at the infinity. 
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Lemma 5.4. The function u(X) — u(X, £) from Proposition \5.!\ exhibits the fol 
lowing behavior as £ — > oo; 



«(A,o = r 1/MW/K ( E%( A )r 2i + o(r 2 "- 2 ) | , 



(5.7) 



where a,j(\) are continuous functions of X and the constant in 0(£ 2 " 2 ) is indepen- 
dent of A 6 A. 



Moreover, the expansion (5/7) can be differentiated in the following sense: 
<9 fc u„ , N <9 fc 



rl/CT - WK ^ a . (A)r2( 

J=0 



-l/a-2n-2-k\ 



(5.8) 



where the constant in 0(£ _1 / cr ~ 2 " _2 ~ i; ) is independent of A G A. 

PROOF: T he proof is similar to bootsrapping arguments used in the regularity 
theory. From (5^) and the assumption (2a + < 2 + 2/ a — d we get immediately 



u(A,0 = 0(r 1/CT ) as £^oo. 
We now rewrite the equation u — Tu as 



u(0 = h + ^i)P(0-\ ± ^P(0 

1 — ie 



t -i ■ • r 

1 — 



(5.9) 

pw^VI 2 ^, 

(5.10) 



where 



7i 



l + i5 
1 - ie 



EW~ x \u? a u. 



Using (5J3) and ( 5.1Cl|) , we obtain (for a suitable ao) 

«(0 = r 1/CT - iw/K (ao + o(r 2 )) 



(5.11) 



We now repeat the procedure and use ( ^10|) with (|5lT|) instead of (|J). After 
integrating by parts in integrals of the form 



rj e 



-urf/2 _ _ 1 



1 

drj 



which come up in the calculation, we get easily 

U(0 = ( ao + aiC 2 + 0(r 4 } ) _ 



Repeating this procedure, we get ( |5.7| ) by indu ction . The expansion (|5.8|) can be 
obtained in a similar way, after we differentiate ( 5.10|) . □ 

Lemma 5.5. The mapping A i— > u(X) from Proposition 5Ji is continuously differ- 
entiable as a mapping from A to X$ 

Proof: We define 4> : A x X® — > X$ as </>(A, u) = u - T(A,w) and we denote 
M = {(A,w(A)), A e A} C A x X$. By using similar estimates as those used in the 
proof of Lemma |5.2|, one can easily see that <j) is differentiable with respect to v and 
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that the (partial) derivative D v 4> is continuous as a map from Ax^ into the space 
C{X$, X#) of bounded linear operators on X#. Moreover, a calculation similar to the 
one leading to (|5.4j) gives 



| J - D„0| < 2C7er P (5-12) 

in A x Bp. The refore, since we assume (^fj), D„0 is invertible for all (A, v) € A x B p . 
By Lemma below, cf> is also differentiable with respect to A at each point of M. and 
the partial derivative T)\<j)(\,v) is continuous in M. Now we can conclude with the 
same arguments as in the Implicit Function Theorem to show that u is differentiable 
and 

||(A) = - [D.„0(A, W (A))]- 1 • [D A 0(A, M (A))] . (5.13) 

□ 



Lemma 5J5 also implies the differentiability of the function A i— > f|(A, £i), since 



the partial derivative of u at £i can be expressed in terms of u due to the fact 
that u solves a differential equation. For example, one can take a smooth func- 
tion (f : [£i,oo) — > R with y>(£i) = 1 and y>(£) = 0, for £ > £i + 1 and write 
= — (u(p)" e?£. Then we express u" from the differential equation for the 
profile ( |l.4D and integrate by parts to obtain 



0£ 

where ip, ip are supported on a finite interval and c\ is a suitable constant. □ 

Lemma 5.6. The mapping <f> : A X X$ — ► X$ defined as cj>(\,v) = v — T(X,v) 
is differentiable with respect to A at each point M. — {(A, u(A)),A G A} C A x X$. 
Moreover, the partial derivative ~D\<fi « s continuous in Ai. 



Proof: Using ( |2.3|) we note that 

^(a,6,z) = z- Q logz(l + 0(|zr 1 )) , (5.14) 

as z — ■> oo in the sector — n/2 < argz < it and a belongs to a compact subset of 
{C | ReC > 0} C C. Let us carry out the proof for the partial derivative d<j>/d\i, for 
example. The other partial derivatives can be handled in a similar way. (If fact, for 
Ao,A2 and A4 the proof is much easier.) When e > one sees easily (using the rapid 
decay of W -1 ^) in that case) that 4> is differentiable in Ai and that 

^-(A, = t|^(A, + f |^(A, £, V) Hv)\ 2a v( V ) dr, . (5.15) 



It follows from (|5.14j) that, for any 1? > 0, the first term on the right-hand side 
is continuous as a mapping from A to X^. For any 62 > 0, the integral on the right- 
hand side has the required continuity properties in the region e% < e < ei, due to 
the rapid decay of VK _1 (r?). The problem is to obtain estimates which are un iform 



for < e < ei. A closer inspection of the integral on the right-hand side of ( |5.15| ) 
reveals that the only terms for which the required continuity is not obvious come from 
differentiating the exponential e~ %Kr> / 2 ( 1_le ) that appears in W~ x (r\) This term can 
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be handled with using integration by parts as follows. Recalling that Ai = n, we must 
estimate an integral of the form 



P(X, tj)t] 



d-i 9 



8 k 



v(rj)\ a v(rj) drj . 



(5.16) 



To estimate this integral we write 

d_ 

8k 



T) d 
2 k drj 



and in tegrate by parts. This eliminates the power rj 1 obtained after differentiating in 
(5.16), and the required estimate follows. □ 



To complete the proof of Theorem 2.1 we need to show that the coordinate 7 = Ao 
can be replaced by j3 = u(X, £1), if r\ and p are chosen properly. We have 



du 
dX~o 



(A,6) = -P(A,6) + 



K(X^i,v)^-H 2<T u(\, V )d V . (5.17) 



(Since Ao is complex, we interpret this, with a slight abuse of notation, as an equation 
between real 2x2 matrices. Another possibility would be to interpret it literally and 
do a similar calculation for d/d Xp.) We carry out the differentiation in the integral 
on the right-hand side and use ( 5.13| ), which gives 



du 
dAo" 



(A) = -[DW(A,u(A))] " P. 



(5.18) 



We note from (5.12) that by a suitable choice of the constant p (to be specified below) 
we can achieve that 



||[D^(A, U (A))]-i <2, 



(5.19) 



i n A w h ere th e norm is taken in the space of linear operators on the spa ce X p . Using 
( 5.17 ) - (|5.19|) , we see easily, by a similar calculation as in the proof of (5^), that by 
a suitable choice of p we can achieve that -§f^{X, £i)/P(A, £1) is close to the identity 
uniformly in A. A suitable choice for the constants p and r\ is, for example, as follows 



p = a£ x 



ri 



where a is sufficiently small. The proof of Theorem 2.1 can now be easily completed. 
□ 
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